Newtons Method
Newton-Raphson method or only Newton's method is one of the more popular methods used for solving nonlinear algebraic equations of the form
. Nonlinear algebraic equations contain powers of variable(s) and/or transcendental functions. One of the major achievements in mathematics was the proof (by Niels Henrik Abel (1802-1829) in 1824 that polynomial equations of degree greater than 4 cannot be solved by means of an algebraic formula (that is, in terms of radicals).
The basis for Newton's method is that the actual root is estimated and the zero of the tangent to the function at that point is determined.
If a real root
is to be assumed for the function
, then one may easily compute ![]()
![]()
![]()
If we draw a line tangent to the curve at ![]()
![]()
![]()
![]()
![]()
![]()
then the tangent line intersects the
-axis at a point
which is expected to be closer to the actual root
than the assumed root
From the figure we see that
solving for
gives
If we repeat the procedure at ![]()
![]()
![]()
![]()
![]()
![]()
and let
be the
-intercept of the second tangent line, we find
Repetitive use of the procedure yields a sequence of approximations that we expect converges to the root
Newton-Raphson Method
Suppose
is differentiable and suppose
represents the unknown root of
Let
denote a number that is chosen arbitrarily as a first guess to
then repetitive use or iteration of
yields a sequence
of approximations that we expect converges to the root
Newton-Raphson Method
Note that Newton's method unfortunately does not always converge to a root if
is not sufficiently close to the actual root.
See also NewtonsMethod .
The formula of a fourth-degree equation is quite complicated as Maple shows below.
| > | restart:MathMaple:-ini(): |
| > | f:=x->x^4-x-1; |
| > | solve(f(x)=0,{x}): |
| > | evalf(%); #evaluates to floating points |
The roots of the equation above can of course be approximated numerically using Maple's fsolve .
| > | fsolve(f(x)=0,{x}); #all real values |
| > | fsolve(f(x)=0,{x},complex); #all values |
The solution of
| > | f:=x->x^5-3*x^3+x^2-23*x+19:
f(x)=0; |
can be found using fsolve .
| > | fsolve(f(x)=0,{x},complex); |
| > | fsolve(f(x)=0,{x}); #real solutions |
| > | plot(f(x),x=-3..3,`f(x)`); |
![]() |
The following shows an iteration process which does not converge.
| > | f:=x->piecewise(x>0,sqrt(x),-sqrt(-x)):
'f(x)'=f(x); |
First using Newton
| > | Newton(f,1,10,0.0001); |
In this case Newton's method fails. The animation below indicates the reason.
| > | NewtonPlot(f, 1, 10,0.0001,grey); |
![]() |
Example 1
a) The volume of a body is
. Find
by solving the equation
.
b) Solve
= 0.
Solution
| > | restart:MathMaple:-ini(): |
a)
We use the formula
|
We start start with
The approximate value is
|
||||||||||||||||||||||
b)
|
The figure shows that there is 3 or 4 roots. Newton method gives
The solution is
The solution is |
||||||||||||||||||||||
The solution is |
The solution is |
||||||||||||||||||||||
Using fsolve , we get
|
|
||||||||||||||||||||||
To go from one approximation to the next approximation can be visualized by animation. |
Click on the figure and start the animation. |
||||||||||||||||||||||
|
Let us start with a larger initial value.
Click on the figure and select Next Frame repeatedly. |